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Abstract 

In this paper, incorporating the property of the vacuum negative pressure, namely, the bag 
constant, we presented a new model of the equation of state (EOS) of quark matter at finite 
chemical potential and zero temperature. By comparing our EOS with Fraga et aVs EOS and 
SQM1 model, one find that our EOS is softer than Fraga et aVs EOS and SQM1 model. The 
reason for this difference is analyzed. With these results we investigate the structure of quark star. 
A comparison between our model of quark star and other models is made. The obtained mass of 
quark star is 1.3 ~ 1.66M and the radius is 9.5 ~ lAKm. One can see that our star's compactness 
is smaller than that of other two models. 
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The research of compact objects is always highly concerned by physicists and astrophysi- 



cists 



2J, because only these objects can test the fundamental principles of physics under 
extreme conditions, for example, high temperature and high density circumstance in the core 
of compact objects. After the conception of quark has been proposed by Gell-Mann and 
Zweig [3|, |4| , some authors pointed out that neutron star may be composed of quark matter 



llj . In Ref. [9], Witten conjectured that strange quark matter may be the true ground 



state of strong interaction. Inspired by this point of view, some authors raised the concep- 
tion of quark star. In order to investigate the structure and property of this kind of star, 
people have developed various models. The equation of state (EOS) of quark matter plays a 
crucial role for determining the star's structure at high density (above the nuclear density: 
2.4 * 10 u g/cm 3 ) and high temperature (T ~ 10 MeV). The theoretical foundation of all 
these models is Quantum Chromodynamics (QCD). However, at present it is not possible to 
obtain reliably the exact EOS of quark matter from first principles of QCD. So, people try 



to find approximate methods which incorporate basic features of QCD. 



bag model 12l-ll5l|. NJL model 16 



19], and perturbative QCD model |20l - 



br example, MIT 
13 . These models 



incorporate some basic features of QCD. For example, the MIT bag model provide mech- 
anism of quark confinement, NJL model can describe dynamical chiral symmetry breaking 
of QCD, while perturbative QCD model work well at high energy scale due to asymptotic 
freedom. But all these models have their own weaknesses. In MIT bag model, quarks in the 
bag are considered as a free Fermi gas. Obviously, MIT bag model violates chiral symmetry 
even in the limit of massless quark. So, some authors have presented some modified bag 
models, for example, chiral bag model 24j, |25|], which uses pion-quark coupling mechanism 
to restore chiral symmetry. The NJL model assumes that the interaction between quarks 
is point-like, so this model is not renormalizable, and it cannot incorporate quark confine- 
ment 2(| . Perturbative QCD cannot be applied to low energy region, and cannot describe 
chiral phase transition and chiral symmetry breaking. From these one can see how difficult 
it is to understand the property of the compact matter. But by the endeavor mentioned 
above, people have learned much valuable information about them. Then it deserves to 
develop some phenomenal methods for studying compact matter. In this paper, we propose 
a renormalizable EOS based on path integral formalism of QCD and incorporate bag model, 
then apply it to investigate the structure of quark star, and compare it with the models 
mentioned above. 
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As is well known, if we know the thermodynamic potential of a system, then all thermo- 
dynamic variables can be determined [27]. So, using various methods to get the system's 



thermodynamic potential is very important. In Refs. [28H36|. one developed a series of 
nonperturbative methods to demonstrate the aspects of strong interaction at finite density 
and zero temperature. Especially, by means of path integral formalism and employing the 



quark propagator proposed in Ref. [37|, a new EOS of QCD at zero temperature and finite 



chemical potential has been proposed by some of the same authors in Ref. [38| . The purpose 
of this paper is to incorporate the property of the vacuum negative pressure, namely, the 
bag constant into this EOS of QCD to give a new model for studying the quark star. In 
order to be self-contained, here let us first give a brief introduction to the EOS of QCD 



proposed in Ref. 



38l | . According to Ref. 



the quark number density reads: 
pGu) = -N c N f Z 2 J j^tr [G r (/i)(p)74] , 



(1) 



where G r (fi)(p) is the renormalized quark propagator at finite chemical potential /z, Z 2 is 
the quark wave-function renormalization constant, N c and Nf denote the number of colours 
and of flavors, respectively, and the trace operation is over Dirac indices. The model quark 
propagator chosen in Ref. [38j is 



a(/i)(p) = V(c 2 ,A 2 )E l . 

j=i \ip + rrij ip + m* 



(2) 



where p = (p,p4 + ifi), rrij = cij + ibj are complex mass scales. If one sets fi = in Eq. 
([2]), one will obtain the quark propagator at zero chemical potential, which is proposed by 
the authors of Ref. [37J under the guidance of the solution of the coupled set of DSEs for 
the ghost, gluon and quark propagator in Landau gauge. The propagator of this form has 
np pairs of complex conjugate poles located at a-,- ± ibj. When some bj is set to zero, the 
pair of complex conjugate poles degenerates to a real pole. The residues Tj are real (note 
ihat a similar meromorphic form of the quark propagator was previously proposed in Ref. 



39[, in which the residues in the two additive terms are complex conjugate of each other). 
The values of parameters in Eq. ([2]) are taken from Ref. [33] (with notations: 2CC stands 
for two pairs of complex conjugate poles, 1R1CC stands for one pair of complex conjugate 
poles and one real pole, and 3R stands for three real poles, respectively). Just as was 



shown in Ref. 



37] . this propagator can be regarded as a good starting point for the study 



of low energy hadron physics because it can simultaneously describe the dynamical chiral 



symmetry breaking and quark confinement which are very important for low energy QCD. 



By means of this propagator, the authors in Refs. |40|, |41( calculated the pion mass and 
decay constant and the quark number susceptibility at finite chemical potential and zero 
temperature, respectively. 

Using path integral methods, one can obtain the following well-known result 



dfj, 



(3) 



where P{ii) is the pressure density. Integrating Eq. ([3]) gives 



P(» = PGu)| M= o+ / dfi'pQl,') 
Jo 



P(p)\, =0 -N c N f Z 2 /V/ j^tr [G r (/i')(p)7 



(2tt)< 



2^/ 
3tt 2 



\ 



a* + J a 2 - + ft 



(4) 



/ 



with 

I(H\OLj,Pj) 



dp! ( p' 2 - -S 



3/2 



4/i 2 



(5) 
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Dq. (jl]) , one has made 
38]). In Ref. Q, the 



where one has defined ay + = m 2 , and in obtaining the last line of 
use of the quark propagator Eq. (jSJ) (for more details, please see Ref. 
constant term P(li = 0) (the pressure of the vacuum) in the EOS was omitted. This does 
not mean that the term P(/i)| At= o is unimportant. In fact, it is an important quantity. It 
enters the energy density, which is relevant for integrating the Tolman-Oppenheimer-Volkoff 
equations. At present it is not possible to calculate reliably 7 :> (/i)| Ai= o from first principles 
of QCD. In this paper, analogous to MIT bag model, we reconsider the effect of this term 
and think that there exists negative pressure at zero chemical potential in vacuum which 



manifests the confinement of QCD. Namely, we identify P(/i)| M= o with —B, where B is the 
vacuum bag constant. In addition, here, for simplicity, we do not consider the chemical, 
thermodynamical equilibrium and electric neutrality conditions, and assume that the star 
is made of pure three flavor quark matter and have no crust structure. Meanwhile we adopt 
the 2CC parameters. The energy density then reads 

dP 

e(ri = -P(p)+p. — + B. (6) 

In this work we take B as an phenomenological model input and choose the following 
three values for B: B = (65.5 Mel/) 4 , (92 MeV) 4 , (108 MeV) 4 . Our EOS is plotted in 



Fig. [T] and Fig. [21 and analogous to Refs. 42j, |43j, we compare our results with other 

| I 

two models: Fraga et a/.'s perturbative QCD model and SQM1 model [44j . From Fig. [1] 
it can be seen that the pressures in our model for three different values of B differ very 
little. It is clear that our EOS is softer than Fraga et a/.'s EOS (A = 2/x) and SQM1 model. 
In the SQM1 model, the value of B is (164.34 MeV) 4 because of the energy limit of 939 
MeV at zero pressure. To improve our EOS's softness, we have to adopt relatively small 



value of B and take it as a phenomenological parameter only. While in Ref . [45j , the value 
of B (namely, (145 MeV) 4 < B < (162 MeV) 4 ) is determined by the requirement that 
the two flavor quark state is unstable though three flavor quark state is the true ground 
state. In Fig. [2] it can be seen more clearly that our EOS are relatively soft, Fraga et a/.'s 
EOS are moderate, while the SQM1 model is the stiffest. This can be easily explained: 
our EOS incorporates confinement and dynamical chiral symmetry breaking, so it is more 
appropriate to describe the property of quark matter at low density region. Fraga et a/.'s 
model is applicable to the massless quark gas in the chiral symmetry restoration region. The 
SQM1 model demonstrates the non-interacting relativistic free quark gas in the bag, so its 
EOS is the stiffest. Applying this EOS, we obtain the structure of quark star by integrating 
the Tolman-Oppenheimer-Volkoff (TOV) equation: 

dP(r) G (e + P)(M + Anr 3 P) 



dr r(r - 2GM) 

dM(r) 



(7) 



4vrr 2 e. (8) 
dr 

The calculated results are shown in Fig. [3] and we also compare Fraga et a/.'s and SQMl's 
results with ours in this figure. For our model one can see that the scaling relation M ~ B~ l l 2 
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FIG. 1. The pressure as function of density, relative to the free quark gas pressure Pf r 
N c N f ^/(l2ir 2 ). In the SQM1, B = (164.34 MeVf 
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FIG. 2. Quark pressure-energy relation. 

is approximately consistent with Witten's result [9|, meanwhile one notes that our radius 
of quark star are larger than that of Fraga et al.'s model and that of SQM1 model. Then 
the density of star in our model is less than that in their models. In our results, according 
to three values of B, the maximum mass, radius and central number density are: mass ~ 
1.75 M , radius ~ 14 Km, Nc ~ 3.2 p for B = (65.5 MeVf, mass -1.42 M , radius ~ 
10.6 Km, Nc ~ 5.2 p for B = (92 MeV) A , mass~ 1.3 M , radius ~ 9.2 Km, Nc ~ 6.6 p 
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FIG. 3. Mass as function of Radius. 



for B = (108 MeV) A , respectively. This phenomenon just reflects the fact that our EOS 
is appropriate to the low density region, while Fraga et a/.'s and SQMl's EOS work well 
at the high density area (namely, chiral symmetry restoration phase and asymptotic free 
region). The mass-energy density relations calculated with and without the presence of B 
are illustrated in Fig. HI The central density calculated with the presence of B is obviously 
larger than the corresponding one calculated without the presence of B. This demonstrates 
that the star's compactness is promoted by adding the bag model constant. In this paper, we 
cannot consider the interior structure problem and assume that the star is homogeneous pure 
three flavors quark matter. Then there is not phase transition between the quark matter and 
the hadronic matter. But that is not imaginable. Actual stars do have inner structure (the 
crust structure has observable effects such as 7-burst, glitch and so on). Meanwhile we do 
not consider the chemical, thermodynamical equilibrium and electric neutrality conditions. 
So our model is only a rough approximation for the study of quark star. Here we should note 
that the main purpose of this work is to study the applicability of our EOS and qualitatively 
compare our results with those of other models. In the future work we should consider these 
factors. 

Although many models of quark star have been proposed in the literature, one may still 
ask whether quark star exists, and if it exists, how to distinguish it from neutron star. In 



Refs. 



46N48| . the authors think that there is no quark matter core in the neutron star. 
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FIG. 4. Mass-Energy Density relation calculated with and without the presence of B 



Particularly, in Ozel's article 49], she concludes that it is impossible for quark core to exist 
in the central part of neutron star. Based on her observational data, she pointed out that the 
limit on the mass is Mass > 2.1 ±0.28 M , the limit on the radius is R > 13.8 ± 1.8Km. 
Therefore, those soft EOS will be ruled out. Of course, many people do not agree with her 
arguments (see, for example, Ref. [50]). On the other hand, the authors in Ref. 51] think 
that CygnusX — 3 may be a strange quark star by measuring high energy cosmic ray from 
it. Usov argues that 7-ray bursters can be an observational signal of strange quark star. 
By the pair creation mechanism of Coulomb barrier in Ref. 52], the computed result of 



emission energy coincides with the observational data. In cloudy bag models [24J, [25], the 
pion decay may be the source of 7-ray burst and provides efficient cooling mechanism. The 
Kepler rotation period (P& may be lower than 1ms) is smaller than the neutron star with 
the same mass }53| . These can provide useful signal to distinguish quark star from neutron 
star. Nowadays we cannot draw a deterministic conclusion as to whether quark star exists, 
because we still do not have full knowledge about QCD and phase transition from nuclear 
matter to quark matter. 

To summarize, in this article, we present a new, renormalizable EOS at finite density and 
zero temperature and consider the effect of vacuum negative pressure. The pressures for 
different parameter B in our model differ very little. The comparison between our EOS and 
the EOSs in other models show that our EOS is softer than other two models. The reason 
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for the difference is discussed. Applying this EOS, one obtains the structure of quark star. 
The results are mass ~ 1.75 M , radius ~ 14 Km, Nc ~ 3.2 p for B = (65.5 MeV) 4 , mass 
~ 1.42A/q, radius ~ 10.6 Km, Nc ~ 5.2 po 

for B = (92 MeV) 4 , mass ~ 1.3 M , radius 
~ 9.2 Km, Nc ~ 6.6 po for S = (108 MeV) A , respectively. By incorporating bag model, 
one can see that the compactness of our previous EOS is improved. 
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